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_ This paper describes a series of numerical experiments with 
ntegration techniques for differential equations with an aim 
develop criteria for choosing numerical integration methods 
‘or astronautical trajectory and orbit computations. 
~ The requirements of these astronautical integrations impose 
extreme restrictions on the integration techniques, from the 
point of view of accuracy and the time required for the inte- 
gration. Since accuracy in numerical integration processes is 
usually obtained only at the expense of speed of integration, 
and speed of integration obtained usually at the expense of 
accuracy, this paper discusses integration methods that enable 
one to make the best of this inevitably difficult situation. 


Introduction 


This discussion deals with various methods of numeri- 
eal integration of differential equations with particular 
concentration on second order differential equations in 
which the first derivatives are absent (as in the case of 
drag-free motion in a gravitational field). The tech- 
niques described herein are presented in standard 
numerical analysis texts; however the numerical com- 
parisons included in this paper are not found in such 
sources. 

Methods that were considered can be placed in three 
veneral categories: 

1) Difference techniques (including sum-difference). 

2) Runge-Kutta and its variations. 

3) Predictor-corrector techniques. 

For the methods involving difference tables, a meas- 
ure of accuracy can be gained by using methods which 
are specifically suited to second order equations with 
first derivatives absent. Such a method is the Second 
Sum (or Gauss-Jackson) Method, which is a variation 
of a direct difference scheme involving only differences of 
the second derivative. This technique has been used 
with great success by astronomers for many years. A 
heuristic discussion of the advantages of this technique 
san be found in reference [1] and its derivation can be 
found in [2]. 

The Runge-Kutta method that was studied was the 
standard fourth order method [3]. The Gill modification 
of Runge-Kutta was not studied since the saving in 
storage space that would result is insignificant in a 
nedium or large size computer, and the gain in accuracy 
vom less round-off error is insignificant, since Runge- 
Kutta should not be used on such a long computation 
hat round-off could become a prominent error source. 
Phe advantages of Runge-Kutta methods make the 
yrocedures very useful for starting an integration and 
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for integrating up to specified positions in the inde- 
pendent or dependent variables, but the computation 
time per step clearly prohibits the use of these methods 
for long computations. 

Predictor-corrector techniques were investigated with 
primary emphasis on a technique which is specifically 
suited to second order differential equations with first 
derivatives absent [4]. 

As will be indicated, this method proved to be the 
most disappointing of those studied since it suffered 
both accuracy and stability troubles when compared 
with the other two categories of integration techniques. 
Recent modifications of these predictor-corrector 
methods have been developed to overcome the in- 
stability problem [6]. However these methods reduce the 
accuracy of the integration (which was not particularly 
good in original form). 


Difference Table Techniques 


Probably the fastest integration methods rely upon a 
table of differences (and possibly sums) of derivative 
values. These techniques do not have the flexibility of 
convenient interval change that is characteristic of 
Runge-Kutta methods, but they do involve significantly 
less computation time to perform a comparable integra- 
tion provided a uniform integration interval is ad- 
visable for the numerical integration. 

The most generally accepted integration methods of 
the difference table type are Adams Method and the 
Adams-Bashforth Method [5]. The second method is 
actually a variation of the first one in which a correction 
formula is used to improve the accuracy (at the expense 
of doubling the computing time). 

The Second Sum Method is compared with both the 
Adams Method and the Adams-Bashforth Method, and 
also with variations on itself in which a correction 
formula is introduced to improve the accuracy of the 
integration (and double the computation time). 

A series of comparisons is made between these 
methods in the next section, as well as comparisons 
with other classes of methods. Unless specifically indi- 
rated the Second Sum Method will not involve the 
correction formula which doubles the computing time. 
In viewing these tables and subsequent conclusions, it 
should be borne in mind that the normal Second Sum 
and Adams Method require about the same computing 
time and the corrected Second Sum Method and Adams- 
Bashforth (which is simply a corrected Adams Method) 
each take about twice as much computing time as 
Adams or the uncorrected Second Sum Method. 
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Prior to comparing these difference table techniques 
with the other methods studied in this report, we will 
first. investigate some comparisons between Adams, 
Adams-Bashforth, and Second Sum Methods. Following 
this, we will consider the question of whether or not a 
corrector should be used with the Second Sum Method. 
These studies will utilize numerical examples of the 
solution of three differential equations which, as will 
be justified in the next section, test for truncation error, 
instability, and a combination of the two effects. 


a F 3 
Solution of = —y (Testing Truncation Only) 
lz 


Analytic Solution Adams Adams-Bashforth =? 
* sin « iy, NI = 2) 4y, At = .2 4y, At = .2 
5 — .95892427 | —.95929783 | — .95892849 | — .95892320 
8 .98935825 .99002624 .98937139 .98935612 
14 .99060736 .99173687 .990651 44 . 99060368 


In this table, 4V represents a numerical solution in- 
cluding all differences up to 4th differences and =? 
indicates the Second Sum Method. The underlined 
figures represent the figures that are in error and give a 
crude but graphic indication of the error. These ab- 
breviations and conventions will be used throughout 
this section and the next section. 

Notice that the Second Sum Method is superior to 
both Adams and Adams-Bashforth Methods despite 
requiring much less time than the Adams-Bashforth 
Method and about the same time as Adams Method. 
Unfortunately, this situation is reversed when instability 
arises. 


: a : Le. 
Solution of = y (Testing Instability only) 
da, 


Adams-Bashforth z2 
4y, At = .2 4y, Al = .2 


Analytic Solu- 


Aiontene Adams 4y, Al = .2 


5 | .67379470 (—2) .67363983 (—2) 
7 | .91188197 (—8) 91305590 (—3) 
11 | .16701701 (—4) 31340425 (—4) 
27 | .18795288 (—11) 


.67379618 (—2) 
. 91188589 (—3) 
. 16701842 (—4) 
.18795768 (—11) 


.67381116 (—2) 
-91311452 (—8) 
84087840 (—4) 


In this example, it is rather clear that the Adams- 
Bashforth computation is far superior to either of the 
other methods, even when the additional computing 
time is taken into account. The fact that the error does 
not grow, even for an extensive computation, indicates 
that the Adams-Bashforth Method is excellent for in- 
tegrations in which instability 1s a problem. 

To study the interaction between instability and 
truncation error, a combination of exponential and 
sinuosoidal functions was used in an attempt to evaluate 
these methods. 

These results are rather disconcerting in the sense 
that none of the methods could successfully integrate the 
third test differential equation for the specified interval 
size. What is significant is that the Second Sum Method 
was no worse than Adams-Bashforth (which was only 
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F 

y d*y 7; ee rays : ; 

Solution of A = —4y (Testing Instability and Truncation 
be 


2 


% dashes ies WRitoen stag nay = Aas vee *S Fe | 
1 380955987 30955477 - 30956019 - 30955989 
3 70259519 (—2) 70301590 (—2) 70261810 (—2) 70264538 (- 
5 |—.64611811 (—2) |—.64631293 (—2) |—.64629283 (—2) |—.64627829 (= 
8 . 33189276 (—3) 36497755 (—3) 36141908 (—3) 85411446 (— 
11.2/—. 13389468 (—4) |—.87215301 (—3) |—.79004810 (—38) — 66102331 G 


slightly better than Adams Method). On this bas 
plus the first test case, it is a rather obvious conclusi 
that the Second Sum Method is excellent for integratu 
in which truncation error due to function oscillation 1s 
dominant error source. 

Since Adams-Bashforth is simply a modification 
Adams Method in which a corrector formula is ust 
and since this corrector greatly improves the stabil 
of the integration method, it would seem that a c 
rector formula should be used with the Second St 
Method. Because the Second Sum Method is weak or 
in its ability to control instability, this would se¢ 
to be a method of obtaining an excellent general pi 
pose integration technique. 

In order to retain the basic advantage of the Seco 
Sum Method in that it does not use the first derivati 
terms (which are absent from the differential equatia 
of the form being studied), the corrector formula: 
based on the fundamental central difference formula: 


2 aS yy a 1 4M 
m= (v6 $+ 8y > Sag ews 
31 6” 289 By ” 
5 
T 69,480 ° 4" — 3,698,800 °% + 


This formula is actually used as the basis of the Seco 
Sum Method and all of the correctors. The second diff 
ence corrector, abbreviated 6°, uses it directly, the fit 
sum corrector, abbreviated S, uses the formula aft 
applying the sum operator once, and the second st 
corrector, abbreviated >, uses the formula after app 
ing the sum operator sunt 

The results obtained from using these correctors 
summarized in the following tables. The same thi 
differential equations used previously were the basis 
comparison. Unfortunately, the correctors are rath 
disappointing in their effect on these trial solutio: 
sufficiently so that the use of a corrector at all is 
questionable value. 

After viewing the table, it is apparent that none 
the correctors are significantly better than the une 
rected form, in fact, all of them are worse than ¢ 
uncorrected form for at least one of the test integ 
tions. 

On this basis, there is no valid argument for usin 
corrector formula with the Second Sum Method. T 
computation time is essentially doubled by the inclusi 
of the corrector without any gain in accuracy. On t 


: Py 
Solution of ore (Testing Truncation Only) 


Analytic polation. sag vase 2, res Concrien a Corrector a bg nC 
— .95892427 — .95892320 — .95892433 — .95892433 — .95892169 

. 98935825 . 98935612 . 98935829 .98935819 . 98935512 

. 99060736 . 99060368 . 99060758 . 99060683 . 99060613 


, Py 
Solution of ae y (Testing Instability Only) 


Analytic Solution =? Uncorrected 
€ 


4yv, At = .2 


D2, =? Corrector 
4v, At = .2 


2, 6? Corrector 


=2, > Corrector 
2 4vV, At = .2 


4v, At=. 


.49787068 (—1) 
.67379470 (—2) 
.91188197 (—3) 


.49798904 (—1) 
67516844 (—2) 
10134813 (—2) 


49787743 (—1) 
67429677 (—2) 
94898098 (—3) 


49787235 (—1) 
67392934 (—2) 
.92185940 (—3) 


49786171 (—1) 
67238650 (—2) 
.80888140 (—3) 


: dty 
Solution of Fe = —4y (Test Instability and Truncation) 
: % Analytic Solution 2, Uncorrected 2, =? Corrector 22, 22 Corrector D2, 62 Corrector 
é* sin x 4y, At = 1 4y, At = .1 4y, At = .1 Avi At nd 
1 . 30955987 . 30955989 .380955985 .380955985 . 30955979 
3 . 70259519 (—2) .702645388 (—2) . 70269610 (—2) . 70270680 (—2) . 70259020 (—2) 
5 — .64611811 (—2) — .64627829 (—2) — .64632691 (—2) — .64639202 (—2) — .64547292 (—2) 
8 .33189276 (—3) 38411446 (—3) .362883831 (—3) 36416646 (—3) .19269703 (—3) 


basis, we must conclude that the Second Sum Method 
should be used without a corrector formula when con- 
ditions indicate the method should be used. 


Comparison of Integration Techniques 


Since the various methods of integration that are 
being compared rely on different analytical bases, it 
would be quite difficult, if not impossible, to compare 
these methods without resort to empirical methods. 
Simply making reference to the stipulated error terms 
is inadequate since these estimates are frequently only 
upper bounds and reveal nothing concerning the 
stability of the method, 1.e. its ability not to magnify 
the errors that have been made previously. 

To compare these methods with reference to their 
ability to follow accurately the solution of a differential 
equation where truncation error is most significant, they 
were all applied to an integration for the sine solution 
of the second order differential equation 

dy 2 

de 
Since all solutions of this differential equation (includ- 
ing error modes) remain bounded, the only important 
source of error in a relatively short integration is 
truncation. The results of this comparison are sum- 
marized in Table I. 

To study these same integration methods with re- 
spect to their stability, they were applied to the differ- 
ential equation 


TABLE I 


(a) 
Solution of = = y (Testing Truncation Only) 
a 


Analytic |e Adams- ‘ Runge- 
x Solution per are eee { Bashforth Dine: Kutta 
sin x cA gaat Ce eimai Ih a. to) ieee Ac = 1 
1 . 84147098 . 84147098 . 84147169 = - 84147120 . 84147100 
5 |—.95892427 |—.95892427 |—.95893843 |—.95892435 |—.95892489 |—. 95892448 
8 - 98935825 - 98935822 . 98938231 - 98935836 . 98935801 - 98935874 
20 - 91294525 . 91294535 . 91299685 | Nonsense _— - 91294640 
Adams- 
Bashforth 
4v, Av = 2 
1 — 
5 — 95892849 
8 . 98937139 
20 = 
2 
dy 
Saees ot 
dx? 


in an attempt to follow its negative exponential so- 
lution. Since this equation also has a positive ex- 
ponential solution, the integration errors can grow ex- 
ponentially if the mtegration technique is potentially 
an unstable one. The results of this study are given in 
Table IL. 

Since most integration problems are not as idealized 
as either of the examples used in developing Table I 
and Table II and involve both oscillating functions 
(which cause large truncation errors) and exponentially 
growing error terms, a third numerical case was studied. 
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TABLE II 


Solution of “ = y (Testing Instability Only) 


| 
f 
| 


fo 


x Analytic Solution é ” |Second Sum 7y, Aw = 1 Adams 4V, Ax = .1 Adorn Bashigrc Nv, Milne Ax = .1 Runge-Kutta Ax = 
49787068 (—1 49787100 (—1) | .49787067 (—1) = .49795190 (—1) | No entries fo 
: .67379470 = 67381876 (—2) | .67379443 (—2) | .67379481 (—2) | .67984299 e oe oa 
7 91188197 (—3) | .91865500 (—3) | .91188146 (—3) | .91188218 (—3) | .13588766 (=) Vv 
11 .16701701 (—4) 11849989 (—3) | .16701687 (—4) .16701711 (—4) Nonsense | 
13 . 22603294 (—5) Nonsense 22603273 (—5) | .22603325 (—5) beat iea 
18 . 15229980 (—7) 15228916 (—7) | .15229724 (—7) Bee ca 
34 17139084 (—14) 17136784 (—14) | .14011242 (—14) Se oe 
50 19287498 (—21) 19283669 (—21) Nonsense 4 
Adams-Bashforth 4V, 
Ax = .2 
3 49787068 (—1) 
5 67379618 (—2) 
z 91188589 (—3) 
11 16701842 (—4) 
13 . 22603535 (—5) 
18 15230223 (—7) 
34 .17139651 (—14) 
50 19288463 (—21) 
TABLE III 


Soluti es 
olution of — = 


— 4y (Testing Instability and Truncation) 


5 Dwaly fie Solution ae am ns see (oe os Adee ae Runge Kae 
1.0 3.0955987 (—1) 3.0955990 (—1) 8.0955477 (—1) 3.0956791 (—1) = a 
2.0 1.2306002 (—1) 1.2306010 (—1) 1.2306172 (—1) 1.2307874 (—1) 1.2306014 (—1) = 
3.0 7.0259519 (—3) 7.0257570 (—38) 7.0301590 (—3) 7.02433386 (—3) 7.0258696 (—3) — 
4.0 —1.3861322 (—2)| —1.3862467 (—2)| —1.3859552 (—2)| —1.3867013 (—2)| —1.3861365 (—2)| —1.3803029 (—2 
5.0 —6.4611811 (—3)| —6.4631004 (—3)| —6.4631293 (—3)| —6.4638852 (—3)| —6.4610166 (—3) = 
6.0 —6.9260173 (—4)| —6.89753877 (—4)| —6.9901350 (—4) 6.9452092 (—4)| —6.9238004 (—4) = 
7.0 5.9909419 (—4) 6.2164638 (—4) 5.9505955 (—4) 5.9440127 (—4) 5.9817605 (—4) = 
8.0 3.3189276 (—4) 3.7709248 (—A4) 3.6497755 (—A) 3.3098463 (—4) 3.2761140 (—4) 3.3102082 (—4 
ee —1.33889468 (—5)| —1.0576953 (—3)| —8.7215301 (—5)| —8.3339644 (—6) 9.1812591 (—5)| —1.77735386 (—é 


This case involved the solutions of the differential 
equation 


4 


Q 


_—_ 
pa ee 


: : f : : *s sin | 
which has solutions which are of the form e+* f ap 
COS 


so that both truncation errors and instability problems 
(due to the presence of positive exponential terms) would 
be present. The numerical integration attempted to 
follow the solution e “ sin x so that instability problems 
became evident immediately when the numerical so- 
lution began to grow at a positive exponential rate. 
This case should display the interaction between these 
two major sources of integration error. The results are 
summarized in Table III. 

As can be seen from viewing these tables, neither the 
Milne predictor-corrector method nor Runge-Kutta 
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could favorable compare with the difference tab! 
techniques. In order that the Milne Method could ge 
comparable accuracy to the Second Sum Method, | 
was necessary to use the same integration interval an 
therefore twice the computing time, and even with th 
interval, it was more unstable than the not-too-stab! 
Second Sum Method. On this basis of mediocre ac 
curacy and great instability, the method was con 
pletely discarded from further consideration. Tk 
Runge-Kutta method, on the other hand, when co1 
sidered on a basis of equal computing time with tk 
difference table methods (which means four times tl 
interval used for Adams and Second Sum, and twic 
the interval used for Adams-Bashforth), had simp! 
atrocious accuracy though it displayed a most interes 
ing ability to follow a solution without increasing i 
error. lor this reason, the Runge-Kutta Method is we 
suited to exploratory computation where accuracy is n 
emportant. In such cases, the ease of starting the integr. 


tion as well as terminating it at an explicit point, would 
gertainly be of definite value 


steps in which the integration interval aHodld. vary 
greatly (as in a space trajectory from earth to moon 
‘computed using perturbational methods), the ability to 
alter the integration interval in an arbitrarily specified 
‘manner, rather than simply halving and doubling, may 
‘significantly reduce the number of integration steps 
required. This feature, plus the advantages gained from 
‘a self-starting integration procedure, in some cases will 
indicate that the Runge-Kutta method is the most 
desirable procedure. 


Conclusions 


Thus, we arrive at the following general conclusions 
concerning integration methods to be used for extensive 
integrations in which extreme accuracy is important: 

1. Where the integration is known to be stable, the 
Second Sum Method is the best procedure to 
employ. 

2. Where the integration is known to be unstable, the 

: Adams-Bashforth Method should be used. 

3. If an integration is known to be only slightly un- 
stable, the Adams-Bashforth Method should be 
used with a high number of differences (like seven) ; 
however, if the integration is highly unstable, 
Adams-Bashforth should be used with a relatively 
small number of differences (like four). 


4. If space problems arise within the computer, Adams 
Method would be a good choice of method as it 
would involve the shortest program of any of the 
difference table methods. In this case, no more than 
four differences should be used since higher differ- 
ences make this method prohibitively unstable. 

5. There is definitely room for better numerical integra- 
tion methods, particularly for the integration of 
unstable equations with oscillatory solutions. 
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Estimate of the Specific lonization Caused 
by Heavy Cosmic Ray Primaries 
in Tissue or Water* 


T. Foelschet 


Abstract 


The specific ionization caused by the primary particles of 
earbon, neon, iron, and niobium along their paths in water is 
calculated by use of the formulae of Bohr, Bethe and Bloch, 
for the average stopping power. The maximum specific ioniza- 
tion and shape and also half width of the maximum on the 
end of the path of the particles is estimated according to the 
considerations about capture and loss of electrons made by 
Bohr, Brunings, Knipp, Teller, and Bell. 


*From the Air Force Missile Development Center Tech- 
nical Note AFMDC-TN-59-7 dated April 1959. 

+ University of Frankfurt-Main, presently on the Scientific 
and Engineering Staff of the Office of the Chief Scientist, Air 
Force Missile Development Center, Holloman Air Force 
Base, New Mexico. 


Purpose 


The purpose of the following calculation is to estimate 
the specific ionization caused by heavy cosmic ray 
primaries along their paths in living tissue or in water, 
and, provided they do not terminate in disintegration 
stars, especially at the end of their paths. The main 
objective of these investigations is to secure an estimate 
of biological effects. The absolute value of specific ioni- 
zation at different energies and the shape and the half 
width of the ionization maximum at the end of the 
particle’s path are of special interest. To correlate the 
variation of these characteristics with atomic number 
the primary particles carbon, neon, iron, and niobium 
are treated as examples. 


on 
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Theory for Stopping of Heavy Ions in Matter 


To find the ionization along the path the main 
problem involved is calculating the stopping power, 
or the average energy loss per unit path length (dH/ dx) 
of heavy ions with known velocity v. This energy loss 
is due to electronic and nuclear collisions with the atoms 
of the stopping medium. Electronic collisions are inter- 
actions of the charged particles with the electrons of 
the atoms of the stopping medium, while the term “nu- 
clear collisions” refers to recoil phenomena occurring 
near the end of the particle’s path, but not to disintegra- 
tions. Having the expression for the stopping power as 
a function of the velocity or energy per nucleon, the 
range versus velocity can be found by integration and 
the specific ionization versus path can be found by 
eliminating the parameter v or the velocity of the 
particle. 


The Formulae of Bethe-Bloch and Bohr for the Stopping 
Power [4| 


Since the general theory of passage of particles of any 
atomic weight through absorbers of any atomic weight 
is very complex, only the formulae and data applicable 
to the special case of heavy primaries penetrating water 
are discussed. Because of mathematical difficulties in- 
volved in accurate collision theory, semi-empirical data 
must be used to estimate the stopping power or specific 
ionization, especially at the very end of the particle’s 
path, or in the low energy range. 

The contribution of the electronic collisions, or the 
electronic stopping power, may be shown in a general 
expression obtained by Bethe and Bloch. 

This equation is: 


a dg Ar e* <> N(Zi")"Ze 


2mv" ie a: 
we ey rv (1 +4 pes ima) 


forv > Zod 


where: 

y is the logarithmic derivative of the I function, 

Ry denotes the real part of y, 

Zi''e is the effective charge of incident particle, 

v is the velocity of incident particle, 

N is the number of stopping atoms per cc, 

Z, is the nuclear charge of stopping atoms, 

I is the average excitation potential of the atoms of 
the stopping material, 

m is electron rest mass, and 

Vx 18 Bohr’s electron velocity in the ground state of 
the hydrogen atom (e’/(h/27)) 


For vu >> © formula (1) is reduced to: 


h/2xr 
di ,,40 Whe 5 Qmv™ 
dx a mv As iin Tae 2) 


which is Bethe’s formula. 
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Wate 2 


since 
h/2a 


liv« 


Li e Ny 
y(1) = —0.577 and Ry (: +7 yh)? ae In 


we have Bohr’s formula 


dE 4ne’ vette 1.123 mv® \ 
~NZ 22-1 
~ de mv? eS “ute Deh 
z 
For relativistic velocities of the eran particle . 
complete Bethe formula is 3 


¢ 


dE _ yee ay 


dx mv? 2, 


2mv" v v 

-{in?™" ~m (1 = 4) 2 
In the case of heavy primaries traversing water the 
Bloch correction of the Bethe formula is significant fol 


the low energy range, because the requirement for the 
Born approximation : 


(2° 


e . 
> Av: = 24, la 
v> 41 VE | Jn (la) 


is not fulfilled at relatively high velocities of the heavy 
particle. For the end of the path Bohr’s classical ap- 
proximation has to be used. The condition (la) is at 
best valid down to energies somewhat above the lowe 


limits #/A corresponding to v = Zw given below: 
Estimated 

ZWwH E/. aie Residual Faued 
nucl in H:O 
aNb% 41 vg 41.3 1000 u 
ogke*6 26 vH 16.6 400 p 
19 N e?0 10 vx 2.47 55 pw 
gC? 6 vA 0.85 15 pw 


Although these residual ranges are less than one pe 
mille of the path of relativistic particles, in this region 
however, the ionization-maximum occurs. 

The Bloch formula (1) can be used for v > Zovzg 
as shown. That means, considering water as absorbe: 
(Z_ oxygen = 8), about down to v = 8 vy correspond 
ing to 1.6 Mev per nucleon. This, for the heavy particle: 
Nb, Fe, Ne, is in the velocity range where the conditio1 
v XK Zyy 1s fulfilled because Zwater < Zi, (and, there 
fore, v = Ze = 8 vn < Zw). Hence below this ve 
locity the classical Bohr approximation is valid, whicl 
holds independent of the condition v > Z.vyz. The ap 
plicability of the Bloch formula for velocities beloy 
v & 1.6 Mev per nucleon is dubious only for the carbo: 
particle, because in this region v > Zev, is not fulfille 


r is the classical approximation (condition v < Zwy & 
5 Mev per nucleon) valid. 

For an estimate of the shape and half width of the 
pecific lonization-maximum and its absolute value, 
“;' furthermore must be known, and, therefore, the 
omplex mechanism of capture and loss of electrons at 


ow particle velocities must be taken into consideration. 
iy 


Japture and Loss of Electrons at the End of Path, Z°' 


as Function of v 


a When passing through matter with velocities high 
mough so that 


v>Z 4D = Zz, 

wr when the particle’s velocity would be greater than 
he orbit velocity of its innermost electrons would be, 
he incident particle is stripped of all its electrons, so 
hat Z{" is equal to the nuclear charge of the particle. 
for decreasing velocity, however, the probability of 
lectron capture by the nucleus increases, whereas the 
robability of electron loss decreases. Therefore, while 
lowing down, the particle fills up its electron shells; in 
hort, Zi’ decreases with decreasing velocity. For 
articles with higher nuclear charge the electron capture 
rocess begins at much higher velocities than with very 
ght particles like protons and alpha particles and, 
herefore, for heavy primaries the decrease of charge 
annot be disregarded for obtaining the stopping power 
s a function of energy or a range energy relationship of 
ny correctness. A first approximation of the quantita- 
ive relation between effective charge and velocity of 
he particle, based on the Fermi-Thomas statistical 
10del of the atom was given by Bohr, who assumed that 
Hl electrons having a lower orbital velocity than the 
article velocity v would be shaken off, and he gave for 
ne effective charge of the particles: 


bove, vx = ¢/(h/2r) is the Bohr velocity of the elec- 
‘on of the H atom in the unexcited state, Z, being the 
tomic number of the particle. Knipp and Teller [2] 
nd Brunings, Knipp and Teller [2a] gave a variation 
f charge with velocity, which is somewhat closer to 
sperimental evidence, setting the critical electron 
elocity v, (all electrons with viorbitay > Ve remaining 
1 the ion) proportional to the particle velocity v, intro- 
ucing the parameter y by the equation 


Un vv 


sing the Fermi-Thomas model of the absorber atom 
rey calculated the effective charge as a function of 
slocity assuming: 

1. That the characteristic velocity v. is equal to the 


root-mean-square velocity of the energetically most 
loosely bound electron or, alternatively, 

2. equal to the root-mean-square velocity of the owter- 

most electron of the particle. 

The parameter y was determined from the best experi- 
mental evidence available, after computing the stopping 
power, as a function of velocity, for different parameters 
y, chosen within certain limits. The variation of effective 
charge with v, was given in graphical form [2a], p. 658. 
With the knowledge of the parameter y this relation- 
ship can be converted into a relation between particle 
velocity and Z{". 

C. I. Bell [8] has published a more detailed theory 
about the capture and loss of electrons of fast ions. He 
obtains values for Z{"' at a certain velocity about 10 
per cent smaller than those of Brunings-Knipp-Teller’s 
second assumption. In order to give order of magnitude 
calculations in closer agreement with the more recent 
theoretical treatments of the subject of charge exchange 
of fast ions in matter, the second assumption of Brun- 
ings, Knipp, and Teller and the Bohr formula [3] for 
calculation of ionization versus energy is used. The re- 
lations between Z{'' and v, which are used are repre- 
sented in Fig. 1. Substituting Z$" into Bohr’s formula, 
[3], we can calculate the energy loss that is due to 
electronic collisions in the low energy region. 


Nuclear Collisions 


The capture and loss of electrons in the slowing down 
process has a marked effect on the mechanism for the 
energy loss of the particle toward the very end of its 
path. In this lowest energy region, collisions between 
the particle and the nuclei of the stopping material are 
more important than the electronic excitation processes. 
An estimate using the formulas and considerations of 
Bohr [4] shows that the contribution of the nuclear 
stopping power measured in per cent of the electronic 
stopping power Is: 


at 0.3 Mev per nucleon: Nb ~ 4%, Fe ~ 2%, Ne, 


C << 0.5 % @ 


at 0.1 es per nucleon: Nb > 100%, Fe — 50%, 
Neu 8 , ON 3% 


We are at 0.3 Mey per nucleon energy in that region 
where the ionization versus path steeply decreases. For 
the rough first estimate of the half width given here it 
is considered as sufficient to calculate the specific ioniza- 
tion versus path down to 0.2 Mey per nucleon, taking 
into consideration only electronic collisions. The given 
numbers show that above this velocity the contribution 
of nuclear collisions to the overall energy loss per unit 
path can be neglected in the first approximation. For 
the residual range below 0.2 Mev per nucleon empirical 
data are used. 
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Approximative Calculation of Specific Ionization 


Along the Path 


If the stopping power dH/dzx, the decrease of the 
kinetic energy per unit path of the primary of the 
atomic weight A, is given as function of the particle’s 
kinetic energy H, the range FR is obtained by: 


a 


/ °  vdv | 4 dh/A 
Jo fdE(v)| aA es (2) 

{ dx \ da \A 
Eliminating the parameter #/A or v from d#/dx and 
R, or plotting d#/dx as a function of the range Rk, we 
obtain the stopping power along the path. Taking into 
consideration that for each ion pair formed in water 
35 ev is consumed, we get the number of ion pairs formed 
by the primary per unit path or the specific ionization 
along the path dividing dH/dx by 35 ev (if dH /dx is 
given in ev per unit length). These 35 ev comprise not 
only the energy for ionization of the stopping atoms, by 
slow and fast secondary electrons, but also the com- 
paratively large amount of energy, which is lost because 
of mere excitation and dissociation of the molecules or 
atoms. 

The basic curves of energy loss dH /dx or specific 
ionization as function of v or H/A (Fig. 2) for the heavy 
ions of Nb, Fe, Ne, C, are constructed from two parts. 
Above v & Z e*/(h/2r), dE/dx is given by the Bethe 
formulae [2] and [2’]. For such high energies the ratio 
of the stopping powers of two particles having different 
nuclear charge Z,, Z;', but travelling with the same 
velocity through the same material is: 


dE/dx awe 


(dB/Jdxy’  (Z,")? 


For high velocities the stopping power of any charged 
particle can be derived from that of the proton in the 
same medium, which is well known. In Fig. 2 the values 
for protons in water given in [5] multiplied with Z,° 
are used. In the same manner, by comparison with 
the travelled distance of a proton, the range difference 
between the high initial velocity v and the lower limit 
Vv, > Zywy of the nucleus Z, can be calculated by the 
integral 


Ape / dé #1) ee ; | ss OH pee 
» @ Vices ( 
proton 


dx 


or equivalently in terms of the kinetic energy per 
nucleon #/A, when A is the mass number of the par- 
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z$'f OF THE PARTICLE AS FUNCTION OF ITS 
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For low energies of the particle the Bohr formula | 
for the stopping power is used, which in terms of ener, 
per nucleon #/A reads for water (Z. = 10, J = 66 & 


Ry r7cff\2 9 i) 3/2 
CS Na (Zi Y tn (ene /A) Nes 
dx E/A pee ecm 


recht . . ° . a 
Z, asa function of v or H/A is inserted according 
the second assumption of Brunings, Knipp and Tel 
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FIG. 3 I : 
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SPECIFIC IONIZATION CAUSED BY THE PARTICLE 


poONG ITS PATH IN WATER. 


(Fig. 1) with the following values for y: 


7 
Nb 0.37 
Fe 0.41 
Ne 0.45 
C 0.60 


For the intermediate energy region in this first ap- 
proximation the Bloch formula for the stopping power 
is interpolated by the dashed lines which fit in the exact 
values at the lower and higher energies (Fig. 2). In 
this way dH/dx is obtained down to energies 0.2 Mev 
per nucleon. The travelled distance, or range difference 


oft cox) a ol dE IAL 
"(# 02 (-# 
dx dx 


between a given initial velocity v and the lower limit 
dy corresponding to 0.2 Mev per nucleon as function of 
the variable velocity v are then calculated by integra- 
tion with an analog computer. To obtain the total range, 
we have to add the empirical residual ranges at 0.2 
Mey per nucleon which amounts, for the different par- 
ticles, to: 


Residual Range at 
0.2 Mev per nucleon 


Nb 20 p 
Fe 2 
Ne 5 pw 
C 4 iu 


For Nb the value for the light fission fragments is used. 
The value of Fe is interpolated because no experimental 
data were available. The values for Nb, Ne, C, are taken 
from Brunings, Knipp and Teller [2a]. From the sum- 
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mation the total range-energy relations listed in Table 
I result. 

Plotting the stopping power or the specific ionization 
of Fig. 2 at a given energy per nucleon versus the total 
range at the same energy per nucleon the specific ioniza- 
tion along the path is obtained (see Fig. 3 and 4). 
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Results 


Figure 3 shows the specific ionization along the total 
path of relativistic particles having initial velocities of 
10 Bev per nucleon. The maximum of the ionization 
at the end of the path is asymmetrical. The ionization 
decreases very steeply at the very end of the path (Fig. 
3 and 4). Even for the heavy primary Nb this steep 
decrease takes place in a residual range of not more 
than 30 pu. In comparison with the width of this de- 
scending part the half width of the maximum is large 

550 » for Nb). The ascent begins slowly even if the 
maximum is a sharp peak. Its half width does not ex- 
ceed the last yu of the particle’s path as is shown in 
Fig. 4. 

To compare the biological effects of heavy primaries 
with the effect of alpha particles, it is of interest to know 
along which part of the primary’s track the specific 
ionization is higher than the maximum specific ioniza- 
tion caused by an alpha particle (max. 10° ion pairs 
per micron, half width 10 uw). Those features of the 
ionization curve, as the absolute values of the maxima, 
its half widths, and those parts of the path where the 
ionization is higher than 10° ion pairs per micron and 
higher than 10° ion pairs per micron are listed in Table 
DG, 

The half width of the maximum increases with 
atomic weight (not uniformly) by the factor 6 if the 
atomic weight increases by the factor 8. The increase 
of the ionization-maximum itself is roughly proportional 
to the atomic weight of the incident particle. 

It must be emphasized that it is the absolute value 
of the maximum which heavily depends on the relation 
between Z}" and v. Using Brunings, Knipp and Teller’s 
first assumption with empirical values of y the calcula- 
tion results in 40 per cent higher maxima. The above 
given values are to be regarded as rough estimates. 

The main purpose in making these calculations is to 
have a survey of the ranges of heavy primaries in tissue 
and to identify that part of the particle’s path along 
which the particles begin heavily to ionize. 

To give an example, in Fig. 3 it is shown that an Fe 
primary having the relatively low velocity correspond- 
ing to 0.5 Bev per nucleon (28 Bey of energy) has a 
range of about 9 cm in tissue and causes a specific 
ionization higher than 0.8 X 10° ion pairs per p along 
the last 2 cm of its track. In comparison with this, a 4 
Mey alpha particle, as it is emitted from naturally radio- 
active substances, has a range of about 30 yu in tissue. 
The effect of the heavy primary in the last 2 em there- 
fore corresponds at least to the effect of 670 such alpha 
particles. 

Actually the specific ionization along this part of the 
track is up to 10 times higher and is distributed over a 
broader volume [6], so that the biological effect. should 
be considerably higher. 
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TABLE I 
Total Average Range in pu (em) 


0.2 0.6 i 2 3 4 5 6 


Nb 120) 125) 1 3065.) 580 95 | 110 | 125 
Fe | 12 | 19 | 24 | 87 | "52 67 | 87 | 104 
Ne Siete 15 a ORD 92 | 125 | 161 
C A |-10..) 18) 417) 7h 4) 10S) 250 196 


Bev 199 50 100 200 500 1000 "}~.2000 


(cm) 
Nb | 405| 1450 | 3925 | 1.41 | 6.42 | 17.9 | 43.7 
Fe | 455| 1860 | 5300 | 2.05 | 9.42 | 26.8 | 65.3 
(cm) 
Ne | 921] 5371 | 1.54 | 5.21 | 23.0 | 65.0 /158 
C | 1440] 8855 | 2.55 | 8.67 | 38.3 |108 263 


TABLE II 4 
Part of the path with ; 
Maximum ionization where a is: 
Element 
Energy . 3 ; 
ver |Rie| twa | 
n 
sc? | 0.5 9 |1.9X 104 9p 85 pu — 
Ne?) 0.7 20» | 3.0 X 104 1952 500 uw — @ 
oeHe®s| 1.5 38 uw 10° 340 pw ileay rere 20 pw 
aNb%| 2.0 70 |1.7 X 105 5504 | >100 cm* | 3754 


* practically. 
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Some Aspects of the Electrical Properties 
of the Upper Atmosphere’ 


: P. A. Goldbergt 


Abstract 


_ Evaluations are presented of electrical properties of the 
D-region of the ionosphere for altitudes of 150,000 to 300,000 
eet. Certain unique regimes of electrical behaviour of this 
s0rtion of the ionosphere are pointed out. Covered are altitude 
wofiles of a.c. electrical conductivity, dielectric constant, 
ndex of refraction, vertical attenuation factor and penetra- 
jon depth, and long oblique path radio signal attenuation for 
requencies from d.c. to 100 megacycles. Applications to space- 
raft environment problems are treated. 


[Introduction 


The altitude of about 65 kilometers marks an im- 
portant transition region for the geophysical system 
and its environment. Below this altitude lies more 
than 99 per cent of the earth’s “atmosphere”. Above 
t is the beginning of the ‘electrical sea”? which sur- 
rounds our planet, reaching an ionization maximum 
ut some 300 km and then extending diffusely through- 
yut extra-terrestrial space. 

The study presented here is a preliminary analytical 
nodel of the electrical properties of the beginning or 
‘surface”’ region of this electrical sea, the D-region of 
he ionosphere. A principal motivation in this study is 
fhe consideration of the electrical environment of 
spacecraft. Such studies previously [1], [2] have con- 
ventrated on the local electrical environment due to 
onization generated in shock waves. This study is 
‘concerned with the properties of the ambient electrical 
vironment due to ionization in the D-region. Other 
ipplications of the work are connected with the prop- 
rties of the terrestrial waveguide such as appear in 
wroblems of Very Low Frequency propagation [3], 
Whistler propagation [4], [5], and Super Low Frequency 
ropagation [6], [7]. 

Figure 1 places in perspective the role of the D-region 
Nn spacecraft environment problems. The figure 
hows the principal ionosphere regions in terms of the 
vertical distribution of electron density. We see that 
il classes of spacecraft can be subjected to the in- 
luence of the D-region. Hypersonic gliders may operate 
ust below, in or above the D-region. High altitude 
esearch rockets and long range missiles pass through it. 
satellites and lunar and planetary probes are above it 


* Presented at the Fifth Annual Meeting of the AAS, 
Vashington, D. C., Dec. 1958. 
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Fie. 1. Typical distribution of free electrons in the upper 
atmosphere. Proposed hypersonic gliders will operate in the 
lower ionosphere. Missiles, rockets, and satellites pass through 
or must “‘see” through this lower, electrically lossy portion of 
the ionosphere. 


for most of their operations. Thus, all such craft 

generally come directly in contact with the D-region 

or must “look” through it. 

The D-region has certain properties quite different 
from the other ionosphere regions. Although the 
region has the lowest electron density it has the highest 
electron collision frequency of all the ionosphere regions. 
This makes for marked electrical dissipative effects in 
certain electrical regimes. In other regimes the region 
has a pronounced ‘metallic’? type of behaviour. A 
specific objective of this study is the pointing up and 
evaluation of such prominent features of the D-region. 

More generally, the objective of this analysis is to 
develop a broad picture of lower ionosphere electrical 
properties for a wide range of frequencies. The analysis 
proceeds along five major steps: 

1. A model of the electronic parameters; the electron 

density and collision frequency, is set up. 

2. The variation with altitude of the conductivity and 
dielectric constant is evaluated for a wide range 
of frequencies. 

. The propagation parameters, the index of refraction 
and the attenuation factor are then evaluated. 

4. The latter results are examined for the two important 

cases of long vertical and horizontal propagation 


OS 


paths in the D-region. 

5. The ranges of magnitude of all these basic properties 
are then summarized to give an over-all electrical 
picture of the D-region. 
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The analysis covers the wide frequency range of 10 
to 10° c/s. To make the problem tractable for the wide 
range of variables thereby encountered, certain limiting 
considerations are applied. The analysis is carried 
out for the case of a day-time, summer, electron density 
model. This is a maximum ionization condition. In 
evaluating the electrical and propagation parameters 
the earth’s magnetic field is neglected. The results are 
nevertheless good for the cases where the electron 
gyro frequency is not important compared to the col- 
lision and signal frequencies. Also, and more generally, 
when the electric vector of the signal is parallel to the 
geomagnetic vector the analysis holds. In the propaga- 
tion effects analyses the signal wavelength is assumed 
small compared to distance changes of local refractive 
index. This ray treatment is valid for the higher fre- 
quencies. For the lower frequencies such as 10 ke 
the analyses still show the trend of effects and provide 
results directly applicable to techniques developed 
for handling such cases, like that by Poeverlein [8]. 


The D-Region Electron Parameters Model 


The two basic electronic factors determining the 
electrical characteristics of the D-region are the number 
of electrons per unit volume, N, and their frequency of 
collision, vy, with air molecules and atoms. These pa- 
rameters have considerable variation in time. The 
D-region ionization is a maximum near mid-day and is 
present only in traces at night. Also, in addition to a 
seasonal variation there is a strong control of electron 
density by solar activity such as flares which produce 
radio blackouts by enhancement of D-region ionization. 

The D-region electrical model developed here is 
based on a “typical” representation of the electron 
density. For the vertical distribution of electron 
density we use the daytime summer noon model of 
Nertney [9]. The electron collision frequency model 
used here is that of Nicolet [10]. Our resultant electro- 
magnetic model will therefore represent maximum 
electrical effects for a geophysically quiet day. The 
electron models used are shown in Fig. 2. 


Ionosphere Electrical Parameters—General 


In developing the picture of the D-region electro- 
magnetic properties we first construct models of the 
conduction and dielectric properties. From these we 
then construct wave propagation models. The formulae 
to be used in these constructions are taken from jono- 
spheric theory. For the background and development 
of such theory the reader is referred to the brief account 
by Stratton [11] or the comprehensive review by 
Mitra [12]. 

Because there is some variation in convention and 
usage of terminology we shall first set down some 
simple definitions. The differences of terminology have 
to do with inter-relations of conduction, convection 
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Fig. 2. Typical D-Region altitude profiles of electronic ga 
parameters. The electron density is relatively small and # 
collision frequency is relatively quite high. Nevertheles 
both parameters can contribute either strongly or weakly ~ 
electrical effects depending on the electrical excitatic 
frequency. 


and displacement currents which lead to variation 
in use of terms such as “complex conductivity” ¢ 
“complex dielectric constant’. 

Our definitions are based on Maxwell’s equation f 
currents 


V X = ta te 4 ta 


where H is the magnetic intensity, and @.a, t-» and 2 
are conduction, convection, and displacement current 
respectively. Then, for an electric field, #, with angul: 
frequency, w, and the ohmic condition, ¢ -~ o(£), 


V X A = cE + Jwece-® 


where o is the conductivity, «, the permittivity ( 
dielectric constant) relative to that of free space, 
(all units in this paper are those of the rationalize 
m.k.s. system). Contributions for all currents whic 
appear in the real number, g, in the ‘in-phase’? membe 
of, will be lumped together and called the “‘electrie 
conductivity”. Contributions to the +90° out-o 
phase member appearing collectively in the real numb. 
e, Will be called the dielectric constant or relative pe 
mittivity. With these definitions, in the expressio1 
below, the conductivity is always positive; the d 
electric constant may be either positive or negatir 
and is always less than unity. 

The local conductivity of an ionospheric region 
given by the expression, 


o = (Ne’v)/mw + v) = (eqw2v)/(@? + 2) | 


where e and m are the electronic charge and mas 
respectively. In the alternative form, w2 = (N e?)/ (me 


is the “characteristic” or “plasma” frequency, — 


: 
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_ Fie. 3. Electrical conductivity of the lower ionosphere for 
Biitudes from 65 to 91 km and for frequencies from DC to 
100 me. Although numerically small, the cumulative effects 
of conductivity can lead to large dieipauen of electrical 
energy up to frequencies of 10 me, and even higher frequencies 
on geophysically active days. 


The dielectric constant is given by the expression 


é = 1 —(Ne?)/(meo)(w? + v) = 1 — w2/(w2 + 2). (2) 


The above expressions assume that effects of any 
ambient magnetic field can be neglected. Remarks on 
magnetic field effects will be deferred until the later 
section on propagation. 


D-Region Electrical Properties 
Electrical Conductivity 


The variation with altitude and frequency of D-region 
eonductivity is shown in Fig. 3. These values were ob- 
tained by applying expression (1) directly to the elec- 
tronic model of Fig. 2. Additional values for frequencies 
higher than the 10° c/s figures shown can be generated 
by simply translating the 10° c/s curves two decades 
to the left for each decade change in frequency. This is 
because at this end of the frequency regime the con- 
ductivity is to excellent approximation given by 
Ne*y/(mw?). At the lowest frequency, o = Ne?/(myv), 
and the structure of the conductivity profile represents 
principally the effect of the falling off of » with altitude. 
(The collision frequency changes by a factor of a 
hundred for only a twentyfold change in electron 
density. ) 

Comparing high and low frequency cases, at the 
highest frequencies, the conductivity oni ~ Nv/o 
and at the lowest frequencies oi ~ N/y, and so 
ri/olo ~ v/w. Therefore the conductivity at high 
frequencies is much lower than that at low frequencies 
and also does not change with attitude much compared 
to the low frequency case. 

A comparison of D-region conductivity with other 
onized gas conditions and conducting substances is 
nade in Table 1. It is seen that the D-region is much 
ess conducting than metals, ionized shock waves, and 
he higher regions of the ionosphere. (The shock wave 
values are estimates based on classical kinetic theory 


TABLE 1 


Comparison of Electrical Conductivities 


D.C. Conductivity in 


Material Goumeo! 


Coppers(202C ieee ere ee een eee 6 X 10’ 
Carboni(2500,C)re, aye LX e10® 
Carbonk(O2@ ire sas: Sane 3 xX 104 
Tellumiuminrins tag ee eee 0X 102 
Shock wave (Mach 20, 56 km)............ 4 xX 10? 
Shock wave (Mach 15, 56 km)............ eer 02 
Shock wave (Mach 13, 56 km)............ 2 X 10! 
Manganese dioxide.....................+- 2X 10! 
Ionosphere, F-region (800 km)............ 2 X< 101 
Dea, Water Sniper Sea Aen ic ye ces 5 X 10° 
Ionosphere, E-region (110 km)............ 2a eam 
Clave > da aes. pat Mn arnt etree >< 10; 
Distilled water (25°C) een es eee ee Ome 
Ionosphere, D-region (90 km).............. ak Se aS 
Ionosphere, D-region (80 km).............. 4X 10° 
DLA Leh a aie. ieee ee wa eet ee LOies 
Ionosphere, D-region (70 km).. Be lOm 
MOVE cae Nae ten ements 5 eee Ade ene a iy <1 
Musediqueartz....cnen nee ee tik nae er eee 2 10m! 


molecular diameters [2].) Nevertheless, as shown later 
although the local conductivity is relatively low the 
cumulative effect for paths through this region are 
quite large for a wide range of field frequencies. In 
addition it will be seen the D-region has a strong metallic 
characteristic at low frequencies in the sense that it 
can have a very high index of refraction and reflec- 
tivity. 


Dielectric Constant 


The dielectric constant model shown in Fig. 4 is 
generated by applying expression (2) to Fig. 2. For 
frequencies above 10° c/s the D-region dielectric 
constant is essentially unity and dielectric effects are 
negligible. This is due to the effect of electron inertia 
at these frequencies and the resulting negligible dis- 
placement currents. At the other extreme, at suffi- 


10° C/S 


ALTITUDE IN KM. 


-4 =3 —3 =I fo) l 
DIELECTRIC CONSTANT, €r 


Fic. 4. Altitude profile of the dielectric constant of the 
lower ionosphere for frequencies from DC to 100 me. 
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TABLE 2 
Comparison of Dielectric Constants 
Dielectric Constant at: 
Material 
102 c/s 10° c/s 108 c/s 

Ionosphere, D-region (90 km)....| —0.88 0.19 1.0 
Ionosphere, D-region (80 km).... 0.98 0.99 1.0 
Ionosphere, D-region (70 km).... 0.99 0.99 1.0 
INGE (ORS, WAU) coon Saco eo weds 1.0 
Aga @-nLSO) abit) kee see rte ila 
Balsa) wood teranamws occ sence 1.4 1.4 1-3 
Waselimeteeanen arcane Gene as ace 2.2 2.2 De 
IRolvstivremenes ava eee cs 2.6 2.6 2.6 
Alf Syal = 7 On eer ORR cnet Pare 4.2 3.5 
IN COPLOM Gy ay See een om 7.5 5.9 3.4 
HB phivlgaleonoler assess: 25 23 
Barnum 0X0 Gees ae metas sree 34 
Waters (SoeO)inentnees ee one 58 57 
Barium titanate. ......20.-:..+- 4400 


ciently low frequencies, the dielectric constant attains 
large negative values. This is again because of electron 
inertia, this time, however, appearing in the form of 
displacement currents lagging the exciting field. 

In Table 2 the dielectric constants of the D-region 
are compared with some other dielectrics. The table 
points up the unique dielectric feature of a plasma of 
e, being less than unity and possibly negative. As has 
been pointed out elsewhere [13] this feature of the 
ionosphere may attain importance in the operation of 
antennas. Because the system radiates into a region of 
modified wave impedance, antenna detuning can occur. 
For the relatively low altitude of the D-region, local 
ionization formed in the hypersonic shock waves of a 
spacecraft can mask D-region ‘‘contact”’ or local effects. 
The ambient D-region local effect can still be important, 
however, for locations on the spacecraft in expansion 
flow and for thin extended antennas about which slip 
or free flow is occurring. 


Ionosphere Propagation Parameters —General 


The effect of the D-region on the propagation of 
electromagnetic fields will be developed below in terms 
of the index of refraction, n, and index of absorption, 
k. These factors are related to ¢ and e, as 


(n — jk)? = 


or explicitly as 


éy — Jao/ (weg) 


2h? = —e, + [e? + (o/weo)?]!” 
= =i = w2/(w? -}- y)| (3) 
FL cog?) (cota?) ots (a/weo)?} 1? 


Bn? = e+ le? + (o/ce0)?}"” 
= [1 = wA/(o? + ¥)] (4) 
+{L = w/(u + PP + (o/ce0)?})”. 
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The full magneto-ionic expressions for the inde 
as developed by Appleton and Hartree [14] cont 
terms involving electron gyrofrequency compone 
transverse and longitudinal to the direction of pro 
gation of a plane wave front. The expressions (3) an 
(4) represent the cases where these gyrofrequené 
terms can be neglected such as for the lower es 
the D-region where the collision frequency effect 
dominant or when the gyro frequency compon aia ; 
small because of wave orientation. 

It should be noted that both k and-n are alwal 
positive although they can become quite small. Thes 
propagation properties and the others evaluate 
in the following sections were obtained by applyin 
expressions (3) and (4) to the electrical models ¢ 
Figs. 3 and 4. 


D-Region Propagation Properties 
Index of Refraction 


The altitude variation of the D-region refractiv 
index for frequencies to 108 c/s is shown in. Fig. ¢ 
We see that in the low frequency regime between I 
to 10? c/s the D-region refractive index becomes ver. 
large, attaining values as great as 400. Because th 
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INDEX OF REFRACTION, n 
Fia. 5. Altitude profile of D-Region index of refraction fe 


frequencies from 10 to 108 ¢/s. At low frequencies such as 1¢ 
c/s, the lower ionosphere is strongly reflecting because of th 
large values and gradients of refractive index. This makes th 
earth-ionosphere system a very good waveguide at sue 
frequencies. 


absorption index (discussed in the next section) 
relatively low, the D-region at these frequencies is thu 
very ‘‘shiny” or highly reflecting. This state of affain 
is similar to free electron gas conditions in metals an 
so the D-region may therefore be appropriately de 
scribed as “metallic” at these frequencies. 

At higher frequencies, the refractive index falls t 
lower values, but does not fall to less than two pe 
cent of unity. It should be noted, however, that 2 
higher altitudes such as in the E and F-regions whet 
the collision frequency becomes very small, the inde 
of refraction can approach (but never becomes le: 
than) the value zero. This latter effect occurs. fc 


FREQUENCY: 


v ATTENUATION, &; IN DB/KM 


Fie. 6. D-Region altitude profile of local signal attenuation 
for frequencies from 10! to 108 c/s. The transition curve of 108 
¢/s corresponds to the approximate condition of maximum 
energy transfer from the applied electromagnetic field to 
neutral gas particles in the ionosphere. Energy transfer is 
most effective where the excitation frequency approximates 
the local electron collision frequency. 


excitation frequencies slightly greater or slightly less 
‘than the plasma frequency, as indicated in expression 


(4). 
Attenuation Factor 


- The local electrical absorption properties of the 
D-region are shown in Fig. 6. For convenience in 
propagation calculations the values are given in terms 
of a’, the attenuation in decibels per kilometer instead 
of absorption index k, where a’ = 8680 kw/c, and c 
is the free space velocity of light. 

At low frequencies like 10* c/s the medium is only 

moderately absorbing. Thus, such waves do not pene- 
trate very deeply into the D-region. Coupled with 
large refractive index and index gradients this causes 
the medium to be an excellent reflector at such fre- 
“quencies. 
_ At frequencies of 10° to 10° ¢/s the excitation fre- 
‘quency is approximately equal to the electron collision 
frequency. Maximum coupling for energy transfer 
‘into conduction heating then occurs and the absorption 
becomes a maximum. 

At frequencies greater than 10° c/s the absorption 
decreases rapidly. Nevertheless, even for, say, the case 
of 107 c/s, where a’ is only about 5 X 10-? db/km, 
certain types of cumulative effects can make absorp- 
tion at high frequencies significant. 


Effects of D-Region on Propagation over 
Long Paths 


Cumulative effects of D-region attenuation over 
long paths are considered in the following two sections. 
Among the problems involved for such long path 
geometries are: communications between spacecraft 
both through and above the D-region; radiation of 
electrical noise from the earth into extra-terrestrial 
space and vice versa; signalling between spacecraft 
and ground stations. A particular problem in regard 
to radiation is that of whistler propagation (1 ke—20 ke) 


via magneto-ionic ducting into extra-terrestrial space. 
The evaluations developed here apply directly to such 
propagation since for the lower D-region the collision 
frequencies are much higher than the magnetic gyro 
frequencies of the electrons. 


VERTICAL PATH ABSORPTIO 
94 
(fh) dh 
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Fig. 7. Total absorption of signals for vertical direction 
through the D-Region for frequencies from 10 to 108 ¢/s. Such 
cumulative effects of energy dissipation can lead to signal 
losses as great as 300 db at 10° c/s. Most natural and artificial 
signals originating on the earth below frequencies of 5 me. 
should be difficult to observe on satellites because of this 
filtering action of the D-Region. Special exceptions are 
ionospherie leakages such as by the proposed whistler ducts. 


Vertical Path Absorption 


The total absorption experienced by signals passing 
vertically through the D-region is shown in Fig. 7. 
The values shown are the integrated values of the 
attenuation factor: 


94 
afro =f al(f,h) dh. 

62 
The values for the low frequencies are shown dotted to 
emphasize that these values represent only absorptive 
losses calculated from a’. There will be, in addition, at 
these frequencies, effects due to loss by scattering 
because of the high reflectivity of the region. At fre- 
quencies greater than 107 c/s the losses become neg- 
ligible (less than one db). 


Horizontal Path Absorption 


Total absorption for horizontal paths 1000 km long 
are shown in ig. 8. The values are calculated on the 
basis of a plane-stratified D-region. This class of path 
is characteristic of spacecraft orientations such as: 
two satellites communicating along a path which nearly 
grazes the earth or is reflected from the earth’s surface 
near grazing incidence. Another orientation in this 
class would be that of two hypersonic gliders ‘‘cruising”’ 
in the lower ionosphere. 

It is seen that now signals at 107 c/s could experience 
losses of the order of 50 db. It is not until frequencies 
as high as 108 c/s are used that losses can fall to less 
than one db (on a geophysically quiet day). 
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Fig. 8. Another example of large cumulative effects due to 
the weakly conducting D-Region. For long horizontal paths 
through the D-Region losses can run greater than 10,000 db 
at frequencies of 100 KC. Low frequency radio astronomy 
and radio probing of the interplanetary medium will be made 
possible by space stations away from the earth’s obscuring 
ionospheric blanket. 


Summary 


Some of the basic electrical properties of the D-region 
have been evaluated with a major emphasis on their 
nature over the wide frequency range of 10 to 108 c/s. 
A number of prominent effects appear which help to 
fill in the picture of significant features of the electro- 
magnetic environment of spacecraft. To sum up 
broadly, these properties and effects are as follows: 


D-Region Properties 


The region occupies the lower 25 km of the ionosphere. 
Electron densities are relatively low and vary only some 
twentyfold. However, the energy dissipating factor, 
the electron collision frequency, is quite large and 
has a one hundredfold range in magnitude. This latter 
characteristic is generally the most dominant factor 
in D-region effects. 

Table 3 gives an orders-of-magnitude summary of 
basic properties of the D-region. The electrical con- 
ductivity is relatively low being only some 10~‘ times 
that of sea water at low frequencies. However, it varies 
through a range of some 10' depending on frequency. 
The dielectric “constant” is fairly constant but de- 
velops negative values at the top of the D-region. 
Although low in value the conductivity of the region is 
the dominant factor in refractive properties at very 
low frequencies, making the region quite shiny then. 


D-Region Effects 


Perhaps the most important effect of the D-region on 
spacecraft electrical environment is in regard to radio 
communications. For long paths through the region 
certain unfavorable orientations of signal wavefronts 
with the geomagnetic field can lead to extremely high 
signal loss. At frequencies of 100 me and greater, 
effects should generally disappear. However, even at 
100 me wave amplitudes may occasionally be reduced 
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Altitude sE,stn kinabeee ser oie een 64-91 
Electron density, N, in e/m’........... 108 — 3 X 10° 


TABLE 3 | 
Range in Electrical Properties of D-Region 
Electron collision frequency, v, in sec!. 4 X 10’ — 4 X 10% 


Decade range (min-max) at frequency of: 


T 


Electrical Property 102 c/s 10° c/s _ er c/s # 


min | max 


Conductivity, o,in (ohm- | 1077 | 10-4 |10-7_| 10~*| 1071 10- 


met)! 
Dielectric constant, «, —1 1 |-—0.1) 1 1 il 
Index of refraction, n 1 | 100 peti) 1 1 


Local attenuation factor, |10-!| 10° |10-! | 10+ | 10-*| 10° 
a’, in db/km 


Horizontal attenuation, as,| 102 | 10% 102 | 104 | 10-* | 10° 
in db/1000 km 
Vertical attenuation, an, 10! 10? 107 4 


in db 


ten and possibly a hundredfold on a geophysicall 
disturbed day. 
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